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Abstract. The usual theory of plasma relaxation, based on the selective decay of 
magnetic energy over the (global) magnetic helicity, predicts a force-free state for 
a plasma. Such a force-free state is inadequate to describe most realistic plasma 
systems occurring in laboratory and space plasmas as it produces a zero pressure 
gradient and cannot couple magnetic fields with flow. A different theory of relax- 
ation has been proposed by many authors, based on a well-known principle of irre- 
versible thermodynamics, the principle of minimum entropy production rate which 
is equivalent to the minimum dissipation rate (MDR) of energy. We demonstrate 
the applicability of minimum dissipative relaxed states to various self-organized 
systems of magnetically confined plasma in the laboratory and in the astrophysical 
context. Such relaxed states are shown to produce a number of basic characteristics 
of laboratory plasma confinement systems and solar arcade structure. 



1. Introduction and Motivation: 

Self-organization is a natural process [Ortolani and Schnack, 1993] in which a con- 
tinuous system evolves toward some preferred states showing a form of order on long 
scales. Examples of such self-organization processes are ubiquitous in nature and in 
systems studied in the laboratory. These ordered states are very often remarkably 
robust, and their detailed structures are mostly independent of the way the system 
is prepared. The final state to which a system evolves generally depends on the 
boundary conditions, inherent geometry of the particular device, but is relatively 
independent of the initial conditions. Self-organization in plasma is generally called 
Plasma Relaxation. The seminal theory for plasma relaxation first proposed by 
Taylor (1974) yields a force-free state with zero pressure gradient. Taylor's theory 
was the first to explain the experimentally observed field reversal of the Reversed 
Field Pinch (RFP). 

Self-organized states with zero pressure gradient are often far from experimen- 
tally realizable scenarios. Any realistic magnetic configuration confining plasma 
must have a non-zero pressure gradient. Extensive numerical simulations [ Sato et 
al., 1996, Zhu et al., 1995] have established the existence of self-organized states of 
plasma with finite pressure, i.e., these states are observed to be governed by the 
magnetohydrodynamic force balance equation. Numerical simulations [Dasgupta et 
al, 1995, Watanabe et al, 1997] and experiments [Ono et al., 1993 ] have established 
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that counterhelicity merging of two spheromaks can produce a Field- Reversed Con- 
figuration (FRC). The PRC has a zero toroidal magnetic field, and the plasma is 
confined entirely by poloidal magnetic field. It has a finite pressure, nonzero per- 
pendicular component of current and is more often characterized by a high value of 
plasma beta. The stability and long life categorize FRC as a relaxed state, which is 
not obtainable from Taylor model of relaxation. Since spheromaks are depicted as a 
Taylor state [Rosenbluth and Bussac, 1979], formation of a FRC from the merging 
of two counterhelicity spheromaks is a unique process from the point of view of 
plasma relaxation, where a non-force free state (hence non- Taylor state) c;nic!rges 
from the coalescence of two Taylor states. In astrophysical context, particularly so- 
lar physics scenario, a 3-D MHD simulation by Amari and Luciano [2000], among 
others, showed that, after the initial lielicdty drive, the final "relaxed state is far 
from the ... linear force-free model that would be predicted by Taylor's conjecture" 
and they suggested to derive an alternative variational approach. Such physical pro- 
cesses call for an alternate model for plasma relaxation to broadly accommodate a 
large number of these observations. 

In this context, an important work by Turner deserves special attention. Turner 
first showed (1986) that a relaxed magnetic field configuration which can sup- 
port/confine a plasma with a finite pressure gradient. Adopting a two-fiuid ap- 
proach, a model of magnetofluid relaxation is constructed for Hall MHD, under 
the assumptions of minimum energy (magnetic plus fluid) with the constraints 
of global magnetic hclicity hybrid helicity, axial magnetic flux and fluid vortic- 
ity. Euler-Lagrangc equations resulted from such variational approach show the 
coupling between the magnetic field and the fluid vorticity. The solutions for such 
coupled equations are shown to be achieved by the linear superposition of the eigen- 
vectors of the 'curl' operator (force- free states). These solutions yield magnetic field 
configurations which can confine plasma with a finite pressure gradient. 

The principle of "minimum rate of entropy production" formulated by Prigogine 
and others [Prigogine, 1947] is believed to play a major role in many problems of 
irreversible thermodynamics. This principle states that for any steady irreversible 
processes, the rate of entropy production is minimum. For most irreversible pro- 
cesses in nature, the minimum rate of entropy production is equivalent to the min- 
imum rate of energy dissipation. A magnetized plasma is such a dissipative system 
and it is appropriate to expect that the principle of minimum dissipation rate of 
energy has a major role to play in a magnetically confined plasma. It is worth men- 
tioning that Rayleigh [1873] first used the term "principle of least dissipation of 
energy" in his works on the propagation of elastic waves in matter. Chandasekhar 
and Woltjer [1958] considered a relaxed state of plasma with minimum Ohmic dis- 
sipation with the constraint of constant magnetic energy and obtained an equation 
for the magnetic field involving a 'double curl' for the magnetic field and remarked 
that the solution of such equation is "much wider" than the usual force free solution. 

Montgomery and Phillips [1988] first used the principle of minimum dissipation 
rate (MDR) of energy in an MHD problem to understand the steady the steady 
profile of a RFP under the constraint of a constant rate of supply and dissipation of 
helicity with the boundary conditions for a conducting wall. Farengo et al [Farengo 
et al, 1994,1995, 2002, Bouzat, 2006] in a series of papers applied the principle of 
MDR to a variety of problems ranging from flux-core spheromak to tokamak. 

In this work, we present an alternative scenario for plasma relaxation which is 
based on the principle of minimum dissipation rate (MDR) rate of energy. This 
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model gives us a non-force free magnetic field, capable of supporting finite pressure 
gradient. We demonstrate the this model can reproduce some of the basic charac- 
teristics of most of the plasma confining devices in the laboratory as well as the 
arcade structure of solar magnetic field. A two-fluid generalization of this model 
for an open system couples magnetic field with flow - so this model is applicable to 
other astrophysical situations. 

The plan of the paper is as follows: In section 2, we present a single fluid MDR 
model for closed system and briefly discuss its applicarions to some of the laboratory 
plasma confinement devices, such as RFP, Spheromak, Tokamak and FRO. We 
present the result of a recent numerical simulation to justify our choice of the energy 
dissipation rate as an effective minimizer in our variation problem. In section 3 we 
describe a generalized version of the MDR model for open two fiuid system and 
its application to solar arcade problem. We summarize a numerical extrapolation 
method for non-force free coronal magnetic field based on our MDR model. Section 
4 concludes our paper. 



2. Single fluid MDR based model 

To obtain both the field reversal and a finite pressure gradient for RFP, Dasgupta 
et al.[ Dasgupta et al., 1998] considered the relaxation of a slightly resistive and 
turbulent plasma using the principle of MDR under the constraint of constant global 
helicity. With the global helicity Km and the dissipation rate R defined as, 

Km= f A-BrfF; R= [ nfdV, (2.1) 

jv Jv 

the variational principle S{R — AK) leads to the following Euler-Lagrange equa- 
tion, 

VxVxVxB = AB. (2.2) 

Solutions of the above equation can be constructed as a superposition of the force 
free equation using Chandrasekhar-Kendall (CK) eigenfunctions [Chandraskhar 
and Kendall, 1957] 

B = aiBi + aaSa + 0:3^3, (2.3) 

where ai,a2,a3 are constants, to be determined from the boundary conditions, 
and 'B'-s are obtained from 

VxBi = XiBi; (i = 1,2,3); A? = A. (2.4) 

We have used the analytic continuation of CK eigenfunctions in the complex do- 
main. However, one can easily ascertain that the resulting magnetic field B is real. 
Moreover, since the superposition of force free fields with different eigenvalues Aj 
is not a force free field, the resultant magnetic field B is not force free, J x B ^ 0. 
An explicit form of the solution in cylindrical geometry with the boundary condi- 
tions for an RFP (perfectly conducting wall) has shown two remarkable features: 
(1) Such a state can support a pressure gradient; and (2) Field reversal is found in 
states that are not force free. 

In Figure 1 a plot of the field reversal parameter F against the pinch parameter 
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Figure 1. Left panel: Plot of F = B^{a)/ (B^) and 9 = Bg{a)/ (B^) (< .. > is the 
volume average and a is the radial distance) showing the field reversal of RFP obtained 
from MDR model. The dotted curve represents the plot for the minimum energy state of 
Taylor. Right Panel: Pressure profile of RFP obtained from MDR model, (/rom Dasgupta 
et al, 1998). 



G is shown. It is found that axial field reverses at a value Q ~ 2.2; and the pressure 
profile of the RFP obtained from our MDR model. 

A success of this theory is its applicability in a tokamak configuration [Bhat- 
tacharya et al., 2000]. In toroidal geometry, the CK eigenfunctions are hypergeo- 
metric functions in the large aspect ratio approximation. Equilibria can be con- 
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structed under the assumption that the total current J = at the edge. The solu- 
tion reproduces the q-profile, the toroidal magnetic field, the pressure profile and 
(7 — J ■'B/B^ for tokamak, which are very close to the observed profiles (e.g., figure 
2). The solutions allow for a tokamak, low-q discharge as well as RFP like behav- 
ior with a change of eigenvalue which would essentially determine the amount of 
volt-scc associated with the discharge characterizing the relaxed state. The poloidal 
Volt-sec/toroidal fltix associated with the discharge can be directly interpreted in 
terms of the parameter Aro {tq is the minor radius of the torus, A = v^) and 
can be shown to increase with Xtq. MDR relaxed state can yield several kinds of 
solutions with distinct q-profiles. For low values of Aro, the q-profile is monotoni- 
cally increasing and hence the solutions resemble a tokamak type discarge together 
with a nearly constant toroidal field, for larger values of Aro the toroidal field re- 
verses at the edge, signiiying an RFP-type behavior and for intermediate values, 
the relaxed states indicate an ultra low-q type of discharge exhibiting a q-profile 
with a pitch minimum and < q < 1. In figure 2, we plot the toroidal magnetic 
field B^,a = J ■ 3/B^,q and pressure profiles for different values of Arc, for a fixed 
value of the aspect ratio = 4.0. The continuous line corresponds corresponds to 
Aro = 1.5, and shows a monotonic increase with r, which is the essential feature of 
tokamak discharge. For this case the the toroidal magnetic field {B^) is essentially 
constant and the pressure profile has a peak at the center of the minor cross section 
of the torus and falls to zero at the edge. A non-constant behavior is shown by 
the a profile. The dotted line corresponds to Aro = 3.0. In this case the toroidal 
field reverses sign at the edge, which corresponds to an RFP like discharge. The 
q-profile shows a similar behavior. The a profile has a dip at r = and exhibits a 
bump near the outer region. Such profiles have been experimentally observed in the 
ETA-BETA-II device [ Antoni et al, 1983]. For the values Arc = 2.5 the q-profile 
shows a nonmonotonic behavior as shown in figure 3. The toroidal magnetic field 
again peaks at the center and slowly decreases towards the edge, showing a para- 
magnetic behavior that is typical of ultra low-q discharge [Yoshida et al., 1986]. 
Pressure profile peaks at the center and goes to zero at the edge. Some of these 
results have been observed in the SINP tokamak [Lahiri et al, 1996] and were also 
reported earlier in REPUTE-1 [Yoshida et al., 1986] and in TORIUT-6 [Yamada 
et al, 1987]. 

This model has been utilized to obtain the spheromak configuration from the 
spherically symmetric solution under appropriate boundary conditions relevant to 
an insulating boundary at the edge [Dasgupta et al., 2002]. These boundary condi- 
tions are are consistent with the classical spheromak solution used by Rosenbluth 
and Bussac [1979]. The most interesting part of this spheromak solution is the non- 
constant profile for J • B/ B^ plotted against normalized distance from the magnetic 
axis. This profile shows peaks outside the magnetic axis, and this feature is quali- 
tatively closer to the experimentally observed profile reported by Hart et al. [Hart 
et al., 1983]. 

The field-reversed configuration (FRC) is a compact toroidal confinement device 
[Tuszewski, 1988] which is formed with primarily poloidal magnetic field and a 
zero or negligibly small toroidal magnetic field. The plasma beta of the FRC is 
one of the highest for any magnetically confined plasma. FRC can be considered 
as a relaxed state [Guo et al, 2005, 2006], but it is a distinctly non- Taylor state. 
The FRC configuration can be obtained for an externally driven dissipative system 
which relaxes with minimum dissipation. In Bhattacharya et al., [2001] we show that 
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Figure 2. q-profile, B^,a, and pressure profile of (i)Tokamak discharge (continuous 
line); (ii) RFP discharge (dotted line); (from Bhattacharya et al., 2000). 
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Figure 3. q-profile, B^^cr, and pressure profile of low-q Discharge; (from Bhattacharya 

et al., 2000). 

the Euler-Lagrangc equations for relaxation under MDR with magnetic energy as 
the constraint represents an PRC configuration for one choice of eigenvalues. This 
solution is characterized by high beta, zero helicity and a completely null toroidal 
field and supports a non-constant pressure profile (non- force free). Another choice 
of eigenvalue leads to a spheromak configuration. A generalization of this work 
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incorporating flow [Bhattacharya et al., 2003] and taking magnetic and flow energies 
and cross helicities as constraints also produces the FRC configuration. This state 
supports field aligned flows with strong shears that may lead to stability and better 
conflnement. 
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Figure 4. Spheromak configuration (left panel) and the non-constant J\\/B derived 
from our model based on MDR (right panel); {from Dasgupta et al., 2002). 
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2.1. Simulation Results: 

To justify the use of the MDR principle in plasma relaxation and the choice of 
the total dissipation rate as a minimizer in the variational principle we performed 
a fully 3-D numerical simulations for a turbulent and slightly dissipative plasma 
using a fully compressible MHD code [Shaikh et al., 2007] with periodic boundary 
conditions. Fluctuations are Fourier expanded and nonlinear terms are evaluated 
in real space. Spectral methods provide accurate representation in Fourier space. 
In our code there is very little numerical dissipation, and hence the ideal invariants 
are preserved. The selective decay processes (in addition to dissipation) depend 
critically on the cascade properties associated with the rugged MHD invariants that 
eventually govern the spectral transfer in the inertial range. This can be elucidated 
as follows [Biscamp, 2003]. Magnetic vector potential in 3D MHD dominates, over 
the magnetic field fluctuations, at the smaller Fourier modes, which in turn leads 
to a domination of the magnetic helicity invariant over the magnetic energy. On the 
other hand, dissipation occurs predominantly at the higher Fourier modes which 
give rise to a rapid damping of the energy dissipative quantity R. A heuristic 
argument for this process can be formulated in the following way. The decay rates 
of helicity K and dissipation rate R = Jy rjj'^dV in the dimensionless form, with the 
magnetic field Fourier decomposed as B(k, t) = Sfebk exp(ik • r), can be expressed 
as. 



— =--E.fcbk, - = - — j:,kW (2.0) 

where S = tr/ta, is the Lundquist number, tr, ta are the resistive and Alfven time 
scales, respectively. The Lundquist number in our simulations varies between 10^ 
and 10^. We find that at scale lengths for which k w S^^"^, the decay rate of energy 
dissipation is ~ 0(1) . But at these scale lengths, helicity dissipation is only ~ 
0(5~^/^) <^ 1. This physical scenario is consistent with our 3D simulations. Figure 
5 shows the time evolution of the global helicity Km, total magnetic energy Wm and 
the total dissipation rate R. It is seen that global helicity remains approximately 
constant, (decaying very slightly during the simulation time) while the magnetic 
energy decays faster than the global helicity, but the energy dissipation rate decays 
even faster than the magnetic energy. The results of the simulation show that total 
dissipation rate can serve as a minimizer during a relaxation process. 



3. MDR based relaxation model for two fluid plasma with 
external drive: 

Recently, Bhattacharya and Janaki [2004] developed a theory of dissipative relaxed 
states in two-fluid plasma with external drive. We describe the plasma (consisting 
of ions and electrons) by the two-fluid equations. The basic advantages of a two- 
fluid formalism over the single-fluid or MHD counterpart are: (i) in a two-fluid 
formalism, one can incorporate the coupling of plasma flow with the magnetic field 
in a natural way, and (ii) the two-fluid description being more general than MHD, 
is expected to capture certain results that are otherwise unattainable from a single- 
fluid description. It should be mentioned that although we arc invoking the principle 
of minimum dissipation rate of energy - instead of the principle of minimum energy 
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Figure 5. Time evolution of decay rates associated with the turbulent relaxation of the 
ideal MHD invariants, viz, global magnetic helicity Km = Jy A ■ BdV, magnetic energy 
W = Jy /2dV and the total dissipation R = Jy riJ^dV, all shown in the same scale. 

- our approach is closely similar to the approach to plasma relaxation including 
Hall terms, as developed by Turner (1986). 

In the two-fluid formalism the canonical momentum Pj canonical vorticity flj, 
for the j-th species are defined respectively as, 



Pj = rrijUj + ^A, ilj = V X Pj, {j = i, e). (3.1) 
A gauge invariant expression for the generalized helicity is introduced by defining 

Kj= [ Pj-^jdr- [ P'j-^'^dV; (3.2) 
Jv Jv 

which is a relative helicity with respect to reference field having P^- as the canonical 
momentum, and Pq, and P^ are different inside the volume of interest and the same 
outside. The gauge invariance condition holds 

{Vlj - ^) • ^ = (B - B') • n = {ujj - uj'j) • n = 0, (cj^ = V x u^), (3.3) 

under the condition that there is no fiow-field coupling at the surface. The above 
equations then represent natural boundary conditions inherent to the problem. 

The helicity balance equation for the generalized helicity, obtained from the bal- 
ance equations for the generalized momentum and vorticity, consists of two terms, 
the helicity injection terms and helicity dissipation term (containing resistivity and 
viscosity). The global helicity can be taken as a bounded function following Mont- 
gomery et al.,(1988) and it is possible to form a variational problem with the helicity 
dissipation rates as constraints. Thus, 

= — 'j>[Ga — Pa • Ua]J^a ' ^^.a — <j> a ' ^a)Ua ' flda 
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(3.4) 



with 




(3.5) 



(3.6) 



where Li = An /.ii / ecnirj is the magnetic Prandtl number for ions, /i^ is the ion 
viscosity and r] is the resistivity. We assumed the electron mass rUe ^ 0. The 
first two terms in the above integral represent the ohmic and viscous dissipation 
rates while the second two integrals represent the ion and electron generalized 
helicity dissipation rates. Aj and Ae are the corresponding Lagrange undetermined 
multipliers. The Euler-Lagrange equations are written as, 



where ip and x s-re arbitrary gauge functions, with 'V'^ip = = 0. Following 
Turner (1986) we can write the exact solution of the coupled equations (3.7) and 
(3.8) as a linear superposition of two CK eigenfunctions, which are the solutions of 
V X Yfc = AfeYfc i.e., 



and a and f3 are constants and quantify the non force- free part, to be obtained from 
some boundary conditions. 

3.1. Application to solar arcade problem 

One of the important features of this formalism is the presence of a nontrivial flow 
and the coupling between the magnetic field and the flow. The solutions of the 
corresponding Euler-Lagrange equation highlight the role of nonzero plasma flow 
in securing a steady MDR relaxed state for an open and externally driven system. 
This model has been very successful in describing the arcade structures observed 
in Solar corona [Bhattacharya et al., 2007]. To model solar arcade- type magnetic 
flelds, we make some simplified assumptions and consider a 2-D magnetic field in 
the x-z plane with the x direction lying on the photosphere surface and z is the 
vertical direction. With appropriate boundary conditions and after some lengthy 
algebra, the solutions for the components of the magnetic fields are obtained as. 




(3.7) 



(3.8) 



B = Yi + aY2 



a;i=Yi+/3Y2 



(3.9) 
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By=0- 



2 cos 7a; — 



cos 



2 cos 7 



: cos(x\/ 7^ — 1) 



(3.10) 
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27 sin 7a; • 



1 cos 7 



sin ■\/7^ — 1 



■ sin(a;-\/7^ — 1) 



with 7 = (fc^ + A^)/A^ (A is the eigenvalue and k is the wave number associated 
with the solution). The same procedures can be adopted for the flow pattern. In 
figure 7 the' magnetic field lines in tlic,r-2 plane is plotted for two different values 
of the parameter 7. The arcade structures are obtained by solving the field line 
equations ^1 = ^ . It is interesting to note that the arcade structure changes with 
the change of the single parameter 7. Figure 6 plots the fiow lines for 7 — 1.5. 
In cylindrical geometry the solutions of eqn. (3.7) can be written in terms of a 
superposition of Bessel functions. Such solutions yield 3-D Flux ropes, which are 
plotted in the Figure (8). 




Figure 8. Plot of 3-D Fluxrope obtained from the solution eqn.(3.7). 



3.2. Non-force Free Coronal Magnetic Field Extrapolation Based on MDR: 

An outstanding problem in solar physics is to derive the coronal magnetic field 
structure from measured photospheric and/or chromospheric magnctograms. A new 
approach to deriving three-dimensional non-force free coronal magnetic field config- 
urations from vector magnetograms, based on the MDR principle has been recently 
developed [Hu, Dasgupta, & Choudhary (2007)]. In contrast to the principle of 
minimum energy, which yields a linear force-free magnetic field, the MDR gives a 
more general non-force free magnetic field with flow. The full MDR-based approach 
requires two layers of vector magnetic field measurements. Its exact solution can be 
expressed as the superposition of two linear force- free fields with distinct A parame- 
ters, and one potential field. The final solution is thus decomposed into three linear 
force- free extrapolations, with bottom boundary conditions derived from the mea- 
sured vector magnetograms, at both photospheric and chromospheric levels. The 
semi-analytic test case shown in Figure 9 illustrates the feasibility and the high 
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performance of this method. This extrapolation is easy to implement and much 
faster than most other nonlinear force-free methods [Hu & Dasgupta 2007] . It also 
takes full advantage of multiple layer solar magnetic field measurements, thus gives 
a more realistic result. 




-0.5 --S.S -0.5 -O.S 

Figure 9. 3D view of magnetic field lines in a box simulating the typical extrapolation 
domain, (a) exact solution; (b) numerical extrapolation result. Bottom image at height 

shows the normal field component distribution with gray scale indicating strongly negative 
(black) to positive (white) field. The great ressemblance between the exact and the numer- 
ical solutions justifies the feasibility and high performance of this method (Hu, Dasgupta, 
& Choudhary 2007). 



4. Conclusion 

We have demonstrated that a relaxation model based on the minimum dissipation of 

energy can successfully explain some of the salient observations both in laboratory 
and space plasma. An MDR based model model in its simpler form is capable of 
yielding basic characteristics of most of the laboratory plasma confinement devices, 
and a generalization of this model for two fluid with external drive can couple 
magnetic field and flow and predict solar arcade structures. The above findings 
definitely points out that the MDR relaxed states applied to astrophysical plasmas 
is a worthy case for further investigations. In future we plan to apply this model 
to investigate magnetic field structure in the sun with an aim to understanding 
coronal mass ejection. 
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1. Introduction 

The layout design for the Journal of Plasma Physics has been implemented as a 
1^1^25 class file. The JPP class is based on the ARTICLE class as discussed in 
the Wr^K. manual (2nd edition). Commands which differ from the standard I^T5t;X 
interface, or which are provided in addition to the standard interface, are explained 
in this guide. This guide is not a substitute for the I^T^jX manual itself. 

Note that the final printed version of papers will use the Monotype Modern 
typeface rather than the Computer Modern available to authors. Also, the measure 
in JPP class is different from Wr^2s article class. For these reasons line lengths 
and page breaks will change and authors should not insert hard breaks in their text. 

1.1. Introduction to BT^ 

The I^T^X document preparation system is a special version of the TJrjX typesetting 
program. I^TJi?^ adds to Tf^ a collection of commands which allow the author to 
concentrate on the logical structure of the document rather than its visual layout. 

I^TJt;X provides a consistent and comprehensive document preparation interface. 
1^1^^ can automatically number equations, figures, tables, and list entries, as well 
as sections and subsections. Using this numbering system, bibliographic citations, 
page references and cross-references to any other numbered entity (e.g. section, 
equation, figure, list entry) are quite straightforward. 

1.2. General style issues 

Use of defaults will result in a pleasing uniformity of layout and font selection. 

Authors should resist the temptation to make ad hoc changes to these. Also avoid 
use of direct formatting unless really necessary. Papers will be edited as usual, and 
this process may be obstructed by the use of inserted line breaks, etc. 

For general style issues, authors are referred to the 'Preparation of manuscripts' 
in the back cover of the journal. The language used in the journal is British English, 
and spelling should conform to this. 

Use should be made of symbolic references (\ref ) in order to protect against late 
changes of order, etc. 

1.3. Submission of 3TgX articles 

Authors who intend to submit a Mjt;X article should obtain a copy of the JPP class 
file. This is available by anonymous FTP from 

ftp . cup . Ccun .ac.uk 

You will find the class file and instructions contained in a single zip archive jpp-cls . zip 
in the directory 

pub/texarchive/ j ournals/latex/ j pp-cls 

The readme . txt (which is the same directory) tells you how to unzip the jpp-cls . zip 
file. There may also be an 'unpacked' directory containing all of the files separately, 
in case of difficulty. If you cannot obtain the JPP files, use the standard article 
class, with the default 'lOpt' option. 

When submitting the final article, ensure that the following are included and are 
clearly labelled. 

{a) A hardcopy printout of the article. 
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(6) The input file (exactly matching the hardcopy), 
(c) A copy of any user-defined macros. 

(rf) If you have used BibT^jX, the .bib, .bbl and .bst files that were used. 

(e) Any other files necessary to prepare the article for typesetting. 

The files for the final article should be supplied on a PC or Macintosh disk. Please 
write on the disk what format it is, as this saves time. Submit the hardcopy and 
disk by post in the normal way. 

1.^. The JPP document class 

The use of document classes allows a simple change of style (or style option) to 
transform the appearance of your document. The CUP JPP class file preserves the 
standard IMi^X interface such that any document which can be produced using the 
standard I^T^jX ARTICLE class can also be produced with the JPP class. However, 
the measure (or width of text) is slightly different from that for ARTICLE; therefore 
line breaks will change and it is possible that equations may need re-setting. 



2. Using the JPP class file 

First, copy the file jpp.cls and upmath.sty (see §2.2) into the correct subdi- 
rectory on your system. The JPP document style is implemented as a complete 
document class, and not a document class option. In order to use the JPP class, 
replace article by jpp in the Xdocumentclass command at the beginning of your 
document: 

\documentclass{article} 

is replaced by 

\documentclass{jpp} 

Some of the standard document class options should be used. Author-defined mac- 
ros should be inserted before \begin{document}, or in a separate file and should 
be included with the submission, see §5. Authors must not change any of the macro 
definitions or parameters in jpp.cls. 

2.1. Document class options 

In general, the following standard document class options should not be used with 
the JPP class file: 

• lOpt, 11 pt and 12pt - unavailable; 

• twoside is the default (oneside is disabled); 

• onecolumn is the default (twocolumn is disabled); 

• titlepage is not required and is disabled; 

• f leqn and leqno should not be used, and are disabled. 

The standard document class options should not be used with the JPP class: 

• lOpt, llpt, 12pt - unavailable. 

• draft, twoside (no associated class file) - twoside is the default. 

• f leqn, leqno, titlepage, twocolumn - unavailable. 
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2.2. Additional 'packages' supplied with jpp . els 

The following additional package ( . sty) files are supplied in the JPP distribution: 

• upmath - provides 'upright' Greek math characters. Requires the AMS amsbsy 

and amsgen packages. 

If your site does not have the AMS Fonts and AMSM^j?^ packages installed, 
we strongly recommend that your site installs them. With them, you can produce 
output which is much closer to the final result. The latest AMS Fonts/ AMS I^T^X 
distributions can be found on your nearest CTAN (Comprehensive TjjpC Archive 
Network) site. 

2.2.1. The UPMATH package The upmath package defines the macros \mathup and 
Xmathbup, which allow access to the symbols in the AMS Euler fonts. 

The upmath package provides macros for upright lower-case Greek (\ualpha- 
\uxi), and for bold lower-case Greek (\ubalpha-\ubxi). The bold upright symbol 
\eta has to be treated differently, in this case use \uboldeta. The upmath package 
also provides \upartial and \ubpartial. 

To use the upmath package, you need to have the AMS eurm/b fonts installed. 

3. Additional facilities 

In addition to all the standard design elements, the JPP class file includes 

the following features: 

• Extended commands for specifying a short version of the title and author(s) 
for the running headlines. 

• Abstract environment 

• Control of enumerated lists. 

• An extended theorem environment, enabling you to also typeset unnumbered 
theorems, etc. 

• Ability to create new theorem-like environments for claims, conjectures, exam- 
ples, problems, remarks, etc. 

• The subequations and subeqnarray environments for sub-numbering equa- 
tions. 

• A proof environment. 

• A acknowledgements environment. 

Once you have used these additional facilities in your document, it can be processed 
only with jpp. els. 

3.1. Titles, authors' names and running headlines 

In the JPP c;lass, the title of the article and the author's name (or aiithors' names) 
are used both at the beginning of the article for the main title and throughout 
the article as running headlines at the top of every page. The title is used on odd- 
numbered pages (rectos) and the author's name (with initials only for first names) 
appears on even-numbered pages (versos). The \pagestyle and \thispagestyle 
commands should not be used. Similarly, the commands \markright and \markboth 
should not be necessary. 

Although the article title can run to several lines of text, the running headline 
must be a single line. Moreover, the title can also incorporate new line commands 
(e.g. \\) but these are not acceptable in a running headline. To enable you to specify 
an alternative short title and author's name, the standard \title and \author 
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commands have been extended to take an optional argument to be used as the 
running headUne: 

\title[A short title] {The full title which can be as long 

as necessary} 

\author [Initials and last neimes of all authors] {The full names of 
all the authors, using letterspacing and \cind before 
the last name in the list} 

\af f iliation{Author ' s affiliation} 

Unlike most other document classes, j pp. els does not follow the convention 
of using the second line of the \author command to write the author's affilia- 
tion. Rather, this must be entered using the \affiliation command as shown 
above. The \and command does not produce separate author /affiliation pairs, but 
only generates the word 'and' in small caps as required by the JPP class. Use the 
predefined macro \ls to get the letterspacing between the letters of the author's 
name. If there is more than one affiliation, use \\ [\af f ilskip] between lines of 
the \af filiation, inserting footnote numbers manually. 

An example will make all this clearer. To produce 

Interesting title 

ALAN N. OTHERS and J.Q. PUBLIC^ 

^University of Moscow, Moscow, Russia 
■^Ngonga University, Nairobi, Kenya 

(Received 12 June 1992) 

you would type 

\title [Short title] {Interesting title} 
\author[A. N. Other and J. Q. Public] 

{A\ls L\ls A\ls N\ns N.\ns 0\ls T\ls H\ls E\ls R$-l$\ns 
\and\ns J.\,Q.\ns P\ls U\ls B\ls L\ls I\ls C$-2$} 
\af f iliation{$"l$University of Moscow, Moscow, RussiaW [\aff ilskip] 

$"2$Ngonga University, Nairobi, Kenya} 
\date{12 June 1992} 

The M^jX \thanks command can be used inside the \author command to insert 
an authors 'Current address' information as a footnote, appearing at the bottom 
of the title page. 

3.2. Abstract 

The JPP class provides for an abstract, produced by 
\begin{abstract} 

\end{abstract} 

This should appear just before the first \section command. 

3.3. Lists 

The JPP style provides the three standard list environments plus an additional 
unmunbered list: 

• Numbered lists, created using the enumerate environment. 
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• Bulleted lists, created using the itemize environment. 

• Labelled lists, created using the description environment. 

The enumerate environment numbers each list item with an italic letter in paren- 
theses; alternative styles can be achieved by inserting a redefinition of the number 
labelling command after the \begin{enumerate>. For example, a list numbered 
with roman numerals inside parentheses can be produced by the following com- 
mands: 

\begin{ enumerate} 

\renewcommaiid{\theenumi}{\romcin{enumi]-} 
\item first item 

\end{ enumerate} 

This produces the following list: 

(i) first item 

(ii) second item 

(iii) etc. 

In the last example, the item alignment is uneven because the standard list in- 
dention is designed to be sufficient for arable numerals rather than the wider roman 
numerals. In order to enable different labels to be used more easily, the eniimerate 
environment in the JPP style can be given an optional argument which (like a 
standard thebibliography environment) specifies the widest label. For example, 

(i) first item 

(ii) second item 

(iii) etc. 

was produced by the following input: 

\begin{enumerate} [(iii)] 
\renewcommEind{\theenumi}{\roman{enumi}} 
\item first item 

\end{enumerate} 

Once you have used the optional argument on the enumerate environment, do 
not process your document with a standard I^TJjX class file. 

3.4- Theorem environments 

The \newtheorem command works as described in the M^TjX manual, but produces 
spacing and caption typefaces required to the JPP style. The preferred numbering 
scheme is for theorems to be numbered within sections, as 1.1, 1.2, 1.3, etc., but 
other numbering schemes are permissible and may be implemented as described in 
the M^jX manual. In order to allow authors maximum flexibility in numbering and 
naming, no theorem-like environments are defined in jpp. els. Rather, you have to 
define each one yourself. Theorem-like environments include Theorem, Definition, 
Lemma, Corollary, and Proposition. For example, 

\newtheorem{lemma}-CLemma} [section] 
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\renewcoiniiiaiid{\theleinma}-{\Romcin{leinma3-} 

The new macro which follows can also produce numbered and unnumbered en- 
vironments. 

3.5. Definitions 

The \newdef inition command can be used to create environments for claims, 
conjectures, examples, problems, remarks, etc. These are typeset in the way as 
theorems, except the text is typeset in roman instead of italic. 

You can have an environment created by \newdef inition number in the same 
sequence as a 'Theorem' by adding the optional [theorem] argument as below: 

\newdefinition{def inition} [theorem] {Definition} 

\begin{def inition} 

This is a definition. 

\end{def inition} 

\begin{def inition*} 

This is an unnumbered definition. 

\end{def inition*} 

which produces the following: 

Definition 3.1. This is a definition. 

Definition. This is an unnumbered definition. 

3. 6. Proof environment 

The standard I^'I]eX constructs do not include a proof environment to follow a 
theorem, lemma etc. (see also §3.4), and so one has been added for the JPP style. 

Note the use in the following examples of an optional argument in square braces 
which may contain any information you may wish to add. For example, 

\begin{theorem} [2] \label{thenv} 

Let the scalar function $T(x,y,t,\biiiath{\omega})$ be a conserved 
density for solutions of {\rm(9)}. Then the two-component function 
\begin-[equation} 

\hbox{\boldit P} = \mathsf -[J}{\cal E} T 
\end{ equation} 

represents the infinitesimal generator of a symmetry group 

for -[\rm (9)}. 

\end{theorem} 

\begin{proof} [Proof of Theorem~\ref {thenv}] 

The assumption about $T$ meams that 

\[ 

\sim \f rac{\partial T}{\partial t} + 
{\cal E}T \bmath{\omega}_t 

= \frac{\partial T}{\partial t} + \{ T, H\}, 
\] 

where $\partial T / \partial t$ refers to explicit 

dependence on $t$. The skew symmetry of $\mathsf{J}$ hence implies 

\begin-[equation} 
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\frac{\partial T}{\partial t> Xsiin \{ H, T\}, 
\end{equation} 

whereupon the operation $\mathsf {J}\cal E$, which commutes with 

$\partial_t$ in its present sense, gives 

\[ 

\frac{\partial \hbox{\boldit P}}{\partial t} 

= \mathsf{J}{\cal E}\{ H, T\}. 

\] 

This equation reproduces the characterisation of symmetries that 
was expressed by (19) , thus showing P to represent a symmetry 
group . 
\end{proof } 

produces the following text: 

Theorem 3.2 (2). Let the scalar function T{x,y,t,u}) be a conserved density 
for solutions of (9). Then the two-component function 

P=]£T (3.1) 

represents the infinitesimal generator of a symmetry group for (9) . 

Proof. [Proof of Theorem 3.2] The assumption about T means that 

dT dT 

where dT/dt refers to explicit dependence on t. The skew symmetry of J hence 
implies 

dT 

-Q^-{H,T}, (3.2) 
whereupon the operation ^£, which commutes with dt in its present sense, gives 

dP 

This equation reproducers the characterisation of symmetries that was expressed by 
(19), thus showing P to represent a symmetry group. □ 

The final □ will not be included if the proof* environment is used. 



4. Mathematics and units 

The JPP class file will insert the correct space above and below displayed maths if 
standard M^rpC commands are used; for example use \ [ ... \] and not $$ ... $$. 
Do not leave blank lines above and below displayed equations unless a new para- 
graph is really intended. 

4-1- Numbering of equations 

The subequations and subeqnarray environments have been incorporated into 

the JPP class file. Using these two environments, you can number your equations 
(1.1a), (1.1&) etc. automatically. For example, you can typeset 

ai = {2nM^/x)^y^ (4.1a) 
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and 



02 = {x/2VL)^ky/M. 



(4.1b) 



by using the subequations environment as follows: 

\begin{subequations} 
\begin{equation} 

a_l \equiv (2\0mega M-2/x) -{\f rac{l}{4}} y-{\f rac{l}{2» 
\end{equation]- 
and 

\begin{equation} 

a_2 \equiv (x/2\0mega) "-CXf rac{l}{2}>k_y/M. 
\end{equation} 
\end{subequations} 

You may also typeset an array such as: 



by using the subeqnarray environment as follows: 

\begin{subeqnarray} 

\dot{Xj k = k \gamma X - \gELinma\delta\eta ,\\ 
\dot{\eta} k = k {\textstyle\f rac{l}-[2}} \delta + 2X\eta . 

\end{ subeqnarray} 

4-2. AMS fonts - especially if you do not have them 

If you need AMS symbols but do not have the AMS fonts you can ensure that 
they will be correctly typeset by taking the following steps. Set up user-defined 
macros that can be redefined by the typesetter to use the correct AMS macros. For 
example, the blackboard bold symbols, sometimes called shell or outline characters, 
are obtained with the AMS macro \mathbb{ . . }. Instead, use a macro definition 
such as: 

y. replace font! 

\newcoinmand\BbbE-[\ensuremath{\mathsf {E}}} % Blackboard bold E 

This substitutes a sans serif character where you want blackboard bold. You can 
typeset the input file and the typesetter is alerted to do the substitution. 

The following example (which uses the Xprovideconmicind macro) will work with- 
out modification by the typesetter, because the Xprovidecommand macro will not 
overwrite any existing \mathbb definition. 

\providecoinincind\mathbb [1] {\ensuremath-C\mathsf {#1}}} 
\newcoinmand\BbbE{\inathbb{E}} Blackboard bold E 

Plain T^TjX provides only \leq and \geq which typeset the Computer Modern 
symbols < and >, respectively. These will be redefined at typesetting to use the 
AMS symbols \geqslant and \leqslant, to give the slanted symbols. 

If you wish to use AMS fonts with WT^^2s you must be using at least version 
2.0. Earlier versions are not supported. 



X = jX — ^St], 



f]=\5 + 2XT^. 



(4.2a) 
(4.2b) 
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4-3. Typeface 

Sometimes, non-italic symbols are required in maths. This section describes how 
these may be obtained using and jpp. els. 

4.3.1. Roman symbols The mathematical operators and constants, such as sin, 
cos, log and cxp (and many others) arc covered by Wr^i^2s macros which ensure 
that they are typeset in roman text, even in math mode: \sin, \cos, \log, \exp. 
Where single letters are concerned (e.g. d, i, e) just use \inathrin in math mode, i.e 
$E=m\mathrm{c}~2$ which typesets as = mc^, giving the correct roman character 
but with math spacing. When the term involves more than one character (e.g. Re 
or Im) text-character spacing is required: 

$\mbox{Re}\;x$ 

which typesets as Re x. 

Where such expressions are used repeatedly, macro definitions can reduce typing 
and editing. The following examples are included in the preamble of the input files 
for this document, jppguide.tex, and of the sample article jppsampl .tex. Authors 

are encouraged to use them and others like them. 

\newcoinincind\Real{\mbox{Re}} °/o do not confuse with TeX's \Re 
\newcoinincind\Iinag{\inbox-[Im}} % do not confuse with TeX's \Iin 

\newcommand\Ai{\mbox{Ai}} 
\newcommand\Bi{\mbox{Bi}} 

4.3.2. Multiletter italic symbols If multiletter symbols are used in maths mode, 
for example Reynolds, Prandtl numbers, etc. the standard maths mode spacing 
between them is too large and text-character spacing is required. As described in 
§4.3.1 (but here for italic letters) use for example 

\newcoininEind\Rey{\mbox{\textit{Re)-}-} 
\newcoinmand\Pran-[\inbox{\textit{Pr}}} 

4.3.3. Upright and sloping Greek symbols Like the italic/roman fonts described 
above, Greek letters for math variables arc printed in the Journal in sloping type, 
whereas constants, operators and units are upright. However, normally pro- 
duces sloping lowercase Greek and upright capitals. The upright lowercase required 
for example for 'pi' and 'mu' (micro) are provided by the upmath package (from 
the AMS Euler fonts). If you don't have the AMS font set, you will have to use the 
normal math italic symbols and the typesetter will substitute the corresponding 
upright characters. 

You will make this much easier if you can use the macros \upi and \umu in your 
text to indicate the need for the upright characters, together with the following 
definitions in the preamble, where they can easily be redefined by the typesetter. 

\providecoinincind\upi{\pi} 
\providecoinmajid\umu{\mu)- 

Notice the use of \providecommand, this stops your definitions overwriting the 
ones used by the typesetter. 

Sloping uppercase Greek characters can be obtained by using the \varGamma etc. 
macros (which are built into the class file). The few required upright Greek capital 
letters (e.g. capital delta for 'difference') are given by the relevant X^T^jX macro, 
e.g. \Delta. 
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4-3.4- Sans serif symbols The \textsf and \mathsf commands change the type- 
face to sans serif, giving upright characters. Occasionally, bold-sloping sans serif is 
needed. You should use the following supplied macros to obtain these fonts. 

\textsf {text} — !■ text \mathsf {math} — > math 

\textsf i{text} — > text \mathsf i{math} — > math 

\textsfb{text} — > text \mathsfb{inath} — > math 

\textsfbi{text} — > text \inathsfbi{inath} — > math 

You can use them like this: 

\newcoinmcind\ssC{\mathsf {C}} 7. for sans serif C 
\newcoinmand\sf sP{\mathsf i{P}} % for sans serif slanted P 
\newcoinmand\sfbsX{\inatlisfbi{X}} % for sans serif bold slanted X 

Note that the bold-slanted font \testsf bi and \mathsf bi use the slanted sans serif 
font cms si because there is no bold-slantcd maths sans serif font in Computer 
Modern! If you use the supplied sans-serif text and math commands the typesetter 
will be able to substitute the fonts automatically. 

4-3.5. Bold math italic symbols If you require bold math italic symbols/letters, 
provides several ways of getting them. Firstly, X^T^jK's \boldmath switch 
which can be used like this: 

$ \mbox{\boldmath $P$} = \mathsf {J}\mathcal{E} T $ 

As you can see it takes quite alot of typing to achieve a bold math italic P. Another 
problem is that you can't use \boldmath in math mode, thus the \mbox forces the 
resulting text into text style (which you may not want). 

You can cut down on the typing by defining the following macro (which is in the 
preamble of this guide): 

\providecommand\boldsymbol [1] {\mbox{\boldmath $#1$}} 

Then for the above example you can define: 

\newcommand\bldP{\boldsymbol{P}} 

the achieve the same result. This still doesn't remove the default to text style 
problem. 

If you have the amsbsy package on your system, then you can remove this lim- 
itation as well by placing a \usepackage{amsbsy} after the \dociimentclass line 

in your document, and then use \boldsymbol for bold math italic symbols (don't 
define \boldsymbol yourself when using cimsbsy). 

4-3-6- Script characters Script characters should be typeset using 2e 's \mathcal 
command. This produces the Computer Modern symbols such as £ and in your 
hard copy but the typesetter will substitute the more florid script characters nor- 
mally seen in the journal. 

4-4- Skewing of accents 

Accents such as hats, over bars and dots are normally centred over letters, but when 

these are italic or sloping greek the accent may need to be moved to the right so that 

it is centred over the top of the sloped letter. For example, \newcommand\hatp{\skew3\hat{p}} 

will produce p. 
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4-5. Units of measure 

Numbers and their units of measure should be typeset with fixed spaces that will 
not break over two lines. This is easily done with user-defined macros. For example, 
52\dynpercm typesets as 52 dynes cm~^, providing the following macro definition 

has been included in the preamble. 

\newcoininEind\dynpercm{\nobreak\mbox{$\ ; $dynes\ , cm$~{-l}$}3- 

5. User-defined macros 

If you define your own macros you must ensure that their names do not conflict 
with any existing macros in plain T^TjX or I^T[t;X (or AMSI^T[t;X if you are using 
this). You should also place them in the preamble to your input file, between the 
\dociimentclass and \begin-C document} commands. 

Apart from scanning the indexes of the relevant manuals, you can check whether 
a macro name is already used in plain TjtjX or I^T^t;X by using the command 
show. For instance, run l^T^X interactively and type \show\<macro_name> at the 
T^pC prompt. (Alternatively, insert the \show command into the preamble of an 
input file and T^TjX it.) 

* \show\Re 

produces the following response from TJn^X: 

> \Re=\mathchar"23C. 
<*> \show\Re 

By contrast \Real is not part of plain T^jX or I^T^jK and \show\Real generates: 

> \Real=undef ined 
<*> \show\Real 

confirming that this name can be assigned to a user-defined macro. 

Alternatively, you can just use Xnewcommand, which will check for the existence 
of the macro you are trying to define. If the macro exists will respond with: 

! LaTeX Error: Command ... already defined. 

In this case you should choose another name, and try again. 

Such macros must be in a place where they can easily be found and modified 
by the journal's editors or typesetter. They must be gathered together in the 
preamble of your input file, or in a separate macros.tex file with the command 
\input -[macros} in the preamble. Macro definitions must not be scattered about 
your document where they are likely to be completely overlooked by the typesetter. 

The same applies to font definitions that are based on Computer Modern fonts. 
These must be changed by the typesetter to use the journal's TJ^X typefaces. In 
this case, you should draw attention to these font definitions on the hard copy that 
you submit for publication and by placing a comment in your input file just before 
the relevant definitions, for example % replace font! 

6. Some guidelines for using standard facilities 

The following notes may help you achieve the best effects with the standard M^iX 
facilities that remain in the JPP style. 
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Table 1. An example table 



Figure 


hA 


hB 




2 


exp (Tlia;) 


exp (niy) 





3 


-1 


exp (Ttia;) 


1 


4 


-4 + 3i 


-4 + 3i 


1.6 


") 


— 2 


— 2 


i.2i 











6.1. Sections 

provides five levels of section headings, only four of which are defined in the 
JPP class file: 

Heading A - \section{. . .} 
Heading B - \subsection{. . .} 
Heading C \subsubsection-[ . . .} 
Heading D \paragraph{. . .} 
There is no subparagraph in the JPP style. 

To obtain non-bold in a bold heading use the usual WTf^2s commands for 
changing typeface; for example \section{Fluctuations in Ca\textsc{ii». 

6.2. Tables 

The table environment is implemented as described in the I^T^riX manual to provide 
consecutively numbered floating inserts for tables. 

JPP class will cope with most table positioning problems and you should not 
normally use the optional positional qualifiers t, b, h on the table environment, 
as this would override these decisions. Table captions should appear at the bot- 
tom of the table; therefore you should place the \caption command before the 
\begin-[tabular>. 

The JPP table environment will insert rules above and below the table as re- 
quired by the JPP style, so you should not attempt to insert these yourself. The 
only time you need to intervene with the rules is when two or more tables fall one 
above another. When this happens, you will get two rules together and too much 
space between the tables; the solution is to add \norule before the end of upper 
table and \f ollowon between tables as shown here: 

\end{center} 

\caption{Dimensionless paraineters}\label-Cdimensionless_p} 

\norule 
\end{table} 

\f ollowon 
\begin{table} 

\begin{center} 

The JPP style dictates that vertical rules should never be used within the body of 
the table, and horizontal rules should be used only to span columns with the same 
headings. Extra space can be inserted to distinguish groups of rows or columns. 

As an example, table 1 is produced using the following commands: 

\begin{table} 

\caption-[An example table} \label-[sample-table} 
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\begin{center} 

\begin-(tabular}-{cccc} 
\hline 

{Figure} & {$hA$} & ■C$hB$} & {$hC$}\\ 
\hline 

2 & $\exp\;(\upi \mathrm{i} x)$ 

& $\exp\; (\iipi \mathrm{i} y)$ & $0$\\ 

3 & $-1$ & $\exp\;(\upi \mathrm{i} x)$ & $1$\\ 

4 & $-4+3\matlirm{i}$ & $-4+3\mathrm{i}$ & 1.6\\ 

5 & $-2$ & $-2$ & $1.2 \mathrm{i}$ 
\end-[tabular}- 

\end{center]- 
\end{table} 

The tabular environment has been modified for the JPP style in the foUowing 
ways: 

(a) Additional vertical space is inserted above and below a horizontal rule produced 
by Xhline 

(6) Tables are centred, and span the full width of the page; that is, they are similar 
to the tables that would be produced by \begin{tabular*}{\textwidth}. 

(c) As with normal I^li^X, if a footnote needs to be inserted into the body of the 
table, the tabular environment should be enclosed in a minipage environment 
of width \textwidth. e.g. 

\begin{minipage}{\textwidth} 
\begin-Ctabular}-[cccc}- 
Xhline 

{Figure} & {$hA$} & {$hB$}\f ootnote{A table must be 
inside a \verb"ininipage" environment if it includes 
table footnotes.} 
& {$hC$}\\ 

\end{tabular} 
\end{minipage} 

Commands to redefine quantities such as \arraystretch should be omitted. If 
the original tabular facilities are needed, there is a new environment 'oldtabular', 
which has none of the reformatting; it should be used in exactly the same way. 

6.3. Illustrations (or figures) 

Artwork for all figures must be supplied separately, to be processed by the Printer, 
and so I^TJt;X's picture environment, for example, cannot be used. An approximate 
amount of space should be left, using the \vspace command. 

The JPP style will cope with most figure positioning problems and you should not 
normally use the optional positional qualifiers t, b, h on the figure environment, 
as this would override these decisions. Figure captions should be below the figure 
itself, therefore the \caption command should appear after the space left for the 
illustration within the figure environment. For example, figure 1 is produced using 
the following commands: 
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Figure 1. An example figure with space for artwork. 

\begin{f igure} 
\vspace{3cm} 

\caption{An example figure with space for artwork. } 
\label{sample-f igure} 
\end{f igure} 

If a figure caption is too long to fit on the same page as its illustration, the caption 

may be typeset as 'Figure X. For caption sec facing page.', and the longer caption 
typeset at the bottom of the facing page. Authors should not concern themselves 
unduly with such details, and may leave pages long. 

6.4- Acknowledgments 

Acknowledgments should appear at the close of your paper, just before the list 
of references and any appendices. Use the acknowledgments environment, which 
simply inserts an 'Acknowledgments' subsection heading. 

There is also an acknowledgment (or acknowledegment) environment for type- 
setting a single acknowledegment. 

6. 5. Appendices 

You should use the standard WI^^ \appendix command to place any Appendices, 
normally, just before the references. This numbers appendices as A, B etc. , equa- 
tions as (Al), (Bl) etc. , and figures and tables as Al, Bl etc. 

If you have only one Appendix, you should use the \oneappendix command 
(instead of \appendix). This command does the same as \appendix, except it 
allows jpp.cls to typeset a single Appendix correctly. 

6. 6. References 

As with standard there are two ways of producing a list of references; either 

by compiling a list (using a thebibliography or theref erences environment), or 

by using BibTjrjX with a suitable bibliographic database. 

References in the text should follow either the Harvard (name/date) or Vancouver 
(numbered) system. In the Harvard system, each reference should be cited in the 
text as 'author(s) (date)' or '(author(s) date)', and the reference list at the end 
of the paper should be in alphabetical order. In the Vancouver system, references 
should be indicated in the text by numbers in square brackets (e.g. [1], [1,2], [1-4], 
etc.), and should be listed at the end of the paper in numerical order of citation. 
For details of the style of entries within the reference list, please see any recent issue 
of the journal. 

6.6.1. References in the text References in the text are given by author and date 
like Dennis (1985) or numbered system [12]. Each entry has a key, which is assigned 
by the author and used to refer to that entry in the text. 
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6.6.2. Harvard (name/date) The following listing shows some references prepared 
in author (date) style of the journal; the code produces the references at the end of 
this guide. 

\begin{thebibliography}{} 

\bibitein[Abramowitz and Stegun (1965)] {AS65} 
{Abramowitz, M. and Stegun, I.\,A.> 1965 
{\ein Handbook of Mathematical Functions}. Dover. 

Xbibitem [Dennis (1985)] {Den85} 
{Dennis, S.\,C.\,R.} 1985 

Compact explicit finite difference approximations to the 
Navier — Stokes equation. In {\em Ninth Intl Conf . on 
Numerical Methods in Fluid DynaimicsX/} (ed. Soubbaramayer 
and J.\,P. Boujot) . Lecture Notes in Physics, vol. 218, 
pp. 23 — 51. Springer. 
\bibitem [Jones (1976)] {Jonl976} 
{Jones, 0.\,C.} 1976 

An improvement in the calculation of turbulent friction in 
rectangular ducts. {\em Trans. ASME\/} J: 
{\em J. \, Fluids Engng\/> {\bf 98}, 173—181. 
\bibitem [Saf f man ( 1990) ] {Saf 1990} 
{Saffman, P.\,G.} 1990 

A model vortex reconnection. {\em J. \, Fluid Mech.} 

{\bf 212>, 395—402. 
\bibitem [Saffman and Schatzman (1982)] {SS1982} 

{Saffman, P.\,G. and Schatzman, J.\,C.} 1982 

Stability of a vortex street of finite vortices. 

{\em J. \, Fluid Mech.} {\bf 117}, 171—185. 
\bibitem [Saffman and Yuen (1980)] {SY1980} 

{Saffman, P.\,G. and Yuen, H.\,C.} 1980 

A new type of three-dimensional deep-water wave of permanent 
form. {\em J. \, Fluid Mech.} {\bf 101}, 797—808. 

\bibitem[Shaqfeh and Koch (1990)] {SK1990} 
{Shaqfeh, E.\,S.\,G. and Koch, D.\,L.} 1990 
Orientational dispersion of fibres in extensional flow. 
{\em Phys. Fluids\/} A {\bf 2}, 1077—1081. 

\bibitem [Wi jngaarden (1968) ] {Wi j 1968} 
{Wijngaarden, L. van} 1968 

On the oscillations near and at resonance in open pipes. 
{\em J.\,Engng Maths\/} {\bf 2}, 225—240 
\bibitem [Williams (1964)] {Wil64} 
{Williams, J.\,A.} 1964 

A nonlinear problem in surface water waves. PhD thesis. 
University of California, Berkeley. 
\end{thebibliography} 

Each entry takes the form 



\bibitem [Author (s) (Date)] {tag} 
Bibliography entry 
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where Author (s) should be the author names as they are cited in the text, Date is 
the date to be cited in the text, and tag is the tag that is to be used as an argument 
for the \cite{} command. Bibliography entry should be the material that is to 
appear in the bibliography, suitably formatted. 

References 

Abramowitz, M. and Stogun, I. A. 1965 Handbook of Mathematical Functions. Dover. 
Dennis, S. C. R. 1985 Compact explicit finite difference approximations to the Navier- 

Stokes equation. In Ninth Intl Conf. on Numerical Methods in Fluid Dynamics (ed. 

Soubbaramayer and J. P. Boujot). Lecture Notes in Physics, vol. 218, pp. 23-51. 

Springer. 

Jones, O. C. 1976 An improvement in the calculation of turbulent friction in rectangular 

ducts. Trans. ASME J: J, Fluids Engng 98, 173-181. 
Saffman, P. G. 1990 A model vortex reconnection. J. Fluid Mech. 212, 395-402. 
Saffman, P. G. and Scliatzman, J. C. 1982 Stability of a vortex street of finite vortices. 

J. Fluid Mech. 117, 171-185. 
Saffman, P. G. and Yuen, H. C. 1980 A new type of three-dimensional deep-water wave of 

permanent form. J. Fluid Mech. 101, 797-808. 
Shaqfeh, E. S. G. and Koch, D. L. 1990 Orientational dispersion of fibres in extensional 

flow. Phys. Fluids A 2, 1077-1081. 
Wijngaarden, L. van 1968 On the oscillations near and at resonance in open pipes. J. Engng 

Maths 2, 225-240 

Williams, J. A. 1964 A nonlinear problem in surface water waves. PhD thesis. University 
of California, Berkeley. 

6.6.3. Vancouver (numbered) system The following listing shows some references 
prepared in Vancouver (numbered) system style of the journal; the code produces 
the references at the end of this guide. 

\begin{theref erences}-[9} 
\bibitem{Jon76} 
{Jones, 0.\,C.}, 

An improvement in the calculation of turbulent friction in 
rectangular ducts. {\em Trans. ASME\/} J: 
{\em J. \, Fluids Engng\/> i\M 98}, 173—181 (1976). 
\bibitem{Saf90} 
{Saffman, P.\,G.>, 

A model vortex reconnection. {\em J. \, Fluid Mech.} 

{\bf 212}, 395—402 (1990). 
\bibitem{SS82} 

{Saffman, P.\,G. and Schatzman, J.\,C.}, 

Stability of a vortex street of finite vortices. 

{\em J. \, Fluid Mech.} {\bf 117}, 171—185 (1982). 
\bibitem{SY80} 

{Saffman, P.\,G. and Yuen, H.\,C.}, 

A new type of three-dimensional deep-water wave of permanent 
form. {\em J. \, Fluid Mech.} {\bf 101}, 797—808 (1980). 
\bibitem{SK90} 

{Shaqfeh, E.\,S.\,G. and Koch, D.\,L.} 

Orientational dispersion of fibres in extensional flow. 
{\em Phys. Fluids\/} A {\bf 2}, 1077—1081 (1990). 
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\bibitem-CWij68} 

{Wijngaarden, L. van}, 

On the oscillations near and at resonance in open pipes. 
{\em JA.Engng Maths\/} {\bf 2}, 225—240 (1968). 
\end{theref erences} 

Each entry takes the form 

Xbibitem-Ctag} 

Bibliography entry 

where tag is the tag that is to be used as an argument for the \cite{} command. 
Bibliography entry should be the material that is to appear in the bibliography, 
suitably formatted. 

References 

[1] Jones, O.C., An improvement in the calculation of turbulent friction in rect- 
angular ducts. Trans. ASME J: J. Fluids Engng 98, 173 181 (1976). 

[2] Saffman, P. G., A model vortex reconnection. J. Fluid Mech. 212, 395-402 
(1990). 

[3] Saffman, P. G. and Schatzman, J. C., Stability of a vortex street of finite vor- 
tices. J. Fluid Mech. 117, 171-185 (1982). 

[4] Saffman, P. G. and Yuen, H. C., A new type of three-dimensional deep-water 
wave of permanent form. J. Fluid Mech. 101, 797-808 (1980). 

[5] Shaqfeh, E. S. G. and Koch, D. L. Orientational dispersion of fibres in exten- 
sional flow. Phys. Fluids A 2, 1077-1081 (1990). 

[6] Wijngaarden, L. van, On the oscillations near and at resonance in open pipes. 
J. Engng Maths 2, 225-240 (1968). 
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Appendix A. Special commands in jpp.cls 

The following is a summary of the new commands, optional arguments and envi- 
ronments which have been added to the standard IM^jX user-interface in creating 
the JPP class file. 



New commands 

\af filiation 
\af f ilskip 

\author 
\f ollowon 
\ls 

\newtheorem 



\newdef inition 



\norule 
\ns 

\oneappendix 

New optional arguments 
[<short title>] 

[<short author>] 



[<widest label>] 

New environments 

acknowl edg ( e ) ment ( s ) 

proof 

proof* 

remark 



subeqnarray 
subequation 
tabular 



use after \author to typeset an affiliation, 
use immediately after \\ in \af filiation to give 
correct vertical spacing between separate affiliations, 
do not use \\ in \author to start an affiliation, 
to remove space between adjacent tables, 
to letter space the author's name, 
this is enhanced so that you can produce unnum- 
bered versions of the environments by using the * 
form. e.g. \begin{theorem*}. 

this environment works like the Xnewtheorem com- 
mand, except it produces environments which are 
typeset in roman instead of italic. Again unnumbered 
'definitions' can be typeset using the * forms, 
to remove the rule below a table, 
to insert space between an author's names, 
does the same as \appendix, except it allows 
jpp. els to typeset a single Appendix correctly. 



in the \title command: to define a right running 
headline that is different from the article title, 
in the \author command: to define a left running 
headline that is different from the authors' names as 
typeset at the article opening. 

in \begin{ enumerate}: to ensure the correct align- 
ment of numbered lists. 



to typeset acknowledgements, 
to typeset mathematical proofs, 
to typeset mathematical proofs without the termi- 
nating proofbox. 

this environment works like the theorem environ- 
ment; it typesets an italic heading and roman text 
to contrast with theorem's small caps heading and 
italic text. 

enables equation numbers in an array to be num- 
bered as (6.1a), (6.2b), etc. 

enables consecutive equations to be numbered (6.1a), 
(6.1b), etc. 

has been modified to insert additional space above 
and below an \hrule and the table caption and body 
is centred with rules full out across the text measure. 
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Appendix B. Notes for editors 

This appendix contains additional information which may be useful to those who are 
involved with the final production stages of an article. Authors, who are generally 
not typesetting the final pages in the journal's typeface (Monotype Modern), do 
not need this information. 

B.0.1. Catchline commands To be placed in the preamble: 

• \date{> 

• \pagerange{} 

• \part{} 

• \pubyear{} 

• \volume{} 

• \doi{> 

B.0.2. Footnotes If a footnote falls at the bottom of a page, it is possible for the 
footnote to appear on the following page (a feature of TJ^X). Check for this. 

B.0.3. Rules between tables If two or more tables fall one above another, add the 
commands \norule and \f ollowon as described in this guide. 

B.0.4- Font substitution Check for use of AMS fonts, bold slanted sans serif, and 

bold math italic and alter preamble definitions to use the appropriate AMS / CUP /Monotype 

fonts for phototypesetter output. 

B.0.5. Font sizes The JPP class file defines all the standard font sizes. For 

example: 

• tiny — This is tiny text. 

• SCriptsize — This is sciiptsizc text. 

• \indexsize — This is indexsize text. 

• \f ootnotesize - This is footnotesize text. 

• \ small - This is small text. 

• \qsmall - This is qsmall text (quotations). 

• \normalsize - This is normalsize text (default). 

• Marge - This is large text. 

• \Large - This is Large text. 

• \LARGE- This is LARGE text. 

All these sizes are summarized in Table B.l. 
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Table B.l. Type sizes for M^jK size-changing commands 



Size 


Size/Baseline 


Usage 


\tiny 


5pt/6pt 




\sc!rip"t siz6 


/ pi / opi 




\indexsize 


8pt/9pt 


catchline. 


\f ootnotesize 


9pt/10pt 


index, footnotes. 


\ small 


9pt/10pt 


'AND' in authors' names, re- 






ceived date, affihation, fig- 






ure and table captions. 


\qsmall 


9.75pt/10.75pt 


quote, quotations. 


Xnormalsize 


10pt/12pt 


main text size, abstract, 




B and C heads. 


Marge 


llpt/13pt 


A head, part no. 


\Large 


14pt/18pt 


part title. 


\LARGE 


17pt/19pt 


article title. 


\liuge 


20pt/25pt 




\Huge 


25pt/30pt 
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An interesting and seminal work on 
various phenomena in fluid mechanics 

ALAN N. OTHERif, 
H.-C. SMITHi and J.Q. PUBLIC^ 

^Department of Chemical Engineering, University of America, Somewhere, IN 12345, 

USA 

^Department of Aerospace and Mechanical Engineering, University of Camford, 
Academic Street, Camford, CF3 5QL, UK 

(Received 1 July 1994 and in revised form 30 October 1998) 

Abstract. Using Stokes flow between eccentric counter-rotating cylinders as a 
prototype for bounded nearly parallel lubrication flow, we investigate the effect of 
a slender recirculation region within the flow field on cross-stream heat or mass 
transport in the important limit of high Peclet number Pe where the enhancement 
over pure conduction heat transfer without recirculation is most pronounced. The 
steady enhancement is estimated with a matched asymptotic expansion to resolve 
the diffusive boundary layers at the separatrices which bound the recirculation 
region. The enhancement over pure conduction is shown to vary as at infinite Pe, 
where is the characteristic width of the recirculation region. The enhancement 
decays from this asymptote as Pe ^ . 



1. Introduction 

The use of integral equations to solve exterior problems in linear acoustics, i.e. to 
solve the Helmholtz equation (V^ -I- = outside a surface S given that 

satisfies certain boundary conditions on S", is very common. A good description is 
provided by 9. Integral equations have also been used to solve the two-dimensional 
Helmholtz equation that arises in water-wave problems where there is a constant 
depth variation. The problem of wave oscillations in arbitrarily shaped harbours 
using such techniques has been examined (see for example Hwang & Tuck 1970; 
Lee 1971). 

In a recent paper 8 have shown how radiation and scattering problems for vertical 
circular cylinders placed on the centreline of a channel of finite water depth can be 
solved efficiently using the multipole method devised originally by Ursell (1950). 
This method was also used by 2 to prove the existence of trapped modes in the 
vicinity of such a cylinder at a discrete wavenumber k < n/2d where 2d is the 
channel width. 

Many water- wave/body interaction problems in which the body is a vertical 
cylinder with constant cross-section can be simplified by factoring out the depth 
dependence. Thus if the boundary conditions are homogeneous we can write the 



f Present address: Plasma Physics Inc., 24 The Street, Lagos, Nigeria 



120 



A.N. Other, H.-C. Smith and J. Q. Public 



velocity potential (j){x,y,z,t) = Re{(j){x,y) cosh k{z + /i)e~"^*}, where the {x,y)- 
plane corresponds to the undisturbed free surface and z is measured vertically 
upwards with z = —h the bottom of the channel. 

Subsequently Callan et al. (1991) proved the existence of, and computed the 
wavenumbers for, the circular cross-section case. It should be noted however that 
experimental evidence for acoustic resonances in the case of the circular cylinder is 
given by Bearman & Graham (1980, pp. 231-232). 

6 provided a theory for determining the trapped-mode frequencies for the thin 
plate, based on a modification of the Wiener-Hopf technique. 

The use of channel Green's functions, developed in § 2, allows the far-field be- 
haviour to be computed in an extremely simple manner, whilst the integral equation 
constructed in § 3 enables the trapped modes to be computed in § 4 and the scatter- 
ing of an incident plane wave to be solved in § 5. Appendix A contains comparisons 
with experiments. 



2. Green's functions 

2.1. Construction of equations 

We are concerned with problems for which the solution, 0, is either symmetric or 
antisymmetric about the centreline of the waveguide, y = 0. The first step is the 
construction of a symmetric and an antisymmetric Green's function, Gs{P, Q) and 
Ga{P,Q)- Thus we require 

(Vl + fc2)G« = (V| + k^)Ga = (2.1) 

in the fluid, where Vp is a gradient operator, except when P = Q, 

Gs,Ga~ l/(27t)lnr asr=|P-Q|^0, (2.2) 

9a 9Ga 



dy dy 
dG, 



on y = d, (2.3) 



= ony = 0, (2.4) 
dy 

Ga on y = 0, (2.5) 

and we require Gg and Ga to behave like outgoing waves as |a;| oo. 
One way of constructing Gs or Ga is to replace (2.1) and (2.2) by 

(V^ + fc2)G, = (V^ + fc2)Ga = -Six - ^)S{y - rj) (2.6) 

and to assume initially that k has a positive imaginary part. 

2.2. Further developments 

Using results from 8 we see that this has the integral representation 
1 f°° 

/ [e-'^'^l^-"' -t- e-'^^^^d-y-r,)! _ Q^y^ ^ ^ (2.7) 

2t Jo 

where 

'^'^ = \ it^-l)K t>l. 



Jq 7 smh k^d 
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In order to satisfy (2.4) wc add to this the function 
1 f °° j^^^^cosh k'^ld — y) 

which satisfies (2.1), (2.3) and obtain 

B{t) = 2e-'=^'* cosh kj{d - ??). (2.8) 

Thus the function 

1 f°° Q-^ld, 

G = -^i{Ho{kr)+Ho{kn))~- j- ———coshkj{d-y)coshk^{d-r,) (2.9) 

7t Jo 7 sinh kja 

satisfies (2.1)-(2.4). By writing this function as a single integral which is even in 7, 
it follows that G is real. 



3. The trapped-mode problem 

The unit normal from D to dD is = {-y'{0),x'{9))/w{9). Now G„ = jYo{kr) + 
Ga where r = {[x{6) — x{ip)]'^ + [y{6) — J/(V')]^}^ and Ga is regular as kr — > 0. In 
order to evaluate dGa{0, 9)/dnq we note that 

as kr — » 0. Expanding x{'4') and y{ijj) about the point ip = 9 then shows that 

^ -[p'{0)p"{e)- p'^{e)-2p'^{e)] asfcr^O. (3.1) 



47tM;3(6')' 
3.1. Computation 

For computational purposes we discretize (3.1) by dividing the interval (0,7t) into 



M segments. Thus we write 

M 



^'^W = ^E'^(^^)a^^«(^'^^X^^)' 0<V<7r, (3.2) 
where 9j = {j — ^)n/M. Collocating at ip = 9i and writing 0^ = 4>{9i) etc. gives 

M 

^•^^ = ^ E i = l,...,M, (3.3) 

where 

dGa{euej)/dn„ i + 3 

\ dG^{eu ei)/dn, + Wl -p\- 2pf]/Anwl i = j. ^ ' ^ 

For a trapped mode, therefore, we require the determinant of the M x M matrix 
whose elements are 

2n 

to be zero. Table 1 shows a comparison of results obtained from this method using 
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Table 1. Values of kd at which trapped modes occur when p{0) = a 



a/d 


M = 4 


M = 8 


Callan et al. 


0.1 


1.56905 


1.56905 


1.56904 


0.3 


1.50484 


1.50484 


1.50484 


0.5 


1.39128 


1.39131 


1.39131 


0.7 


1.32281 


1.32287 


1.32288 


0.9 


1.34479 


1.35160 


1.35185 



Figure 1. Trapped-mode wavenumbers, kd, plotted against a/d for three ellipses: 
, b/a = 0.5; , b/a = 1; , b/a = 1.5. 

two different truncation parameters with accurate values obtained using the method 
of Callan et al. (1991). 

An example of the results that are obtained from our method is given in figure 1. 

Figure 2(a,b) shows shaded contour plots of (j) for these modes, normalized so that 
the maximum value of (j) on the body is 1. 

3.2. Basic properties 
Let 

pi= lim p{x,C}, pu = lim p{x,C) (3.5a,b) 

be the fluid densities immediately below and above the cat's-eyes. Finally let po 
and A^o be the constant values of the density and the vorticity inside the cat's-eyes, 
so that 

{p{x,0,(l>aix,0) = {po,No) ioT Zlix) < C < Z^{x). (3.6) 
The Reynolds number Re is defined by Ut-H/v {y is the kinematic viscosity), the 
length given in wall units is denoted by ( )+, and the Prandtl number Pr is set 
equal to 0.7. In (2.1) and (2.2), and are 

Tij = [uiUj - UiUj) + {uiU^'^^ + uf^^Uj) + uf'^^uf'^^, (3.7a) 

= (^ - Uje) + {UjeSGS + u^GSQ-^ + uSGS0SGS_ (3 Jb) 
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Figure 2. Sliaded contour plots of the potential (f> for the two trapped modes that exist 
for an ellipse with a/d = 1.5, b/d = 0.75. (a) Symmetric about x = 0, kd = 0.96; (6) 
antisymmetric about x = 0, kd = 1.398. 

3.2.1. Calculation of the terms The first terms in the right-hand side of (3.5a) 
and (3.5b) are the Leonard terms exphcitly calculated by applying the Gaussian 
filter in the x- and ^-directions in the Fourier space. 

The interface boundary conditions given by (2.1) and (2.2), which relate the dis- 
placement and stress state of the wall at the mean interface to the disturbance 
quantities of the flow, can also be reformulated in terms of the transformed quanti- 
ties. The transformed boundary conditions are summarized below in a matrix form 
that is convenient for the subsequent development of the theory: 



Qc = 



o w 



iV(a*+u;-iO 



i/3 



Rx V'' 



-(a*a;)-i 






(3i?,-i+c*(ia*)-i) 








-(ia*i?5)-i 















(3.8) 

is termed the displacement-stress vector and Qc the flow-wall coupling matrix. 
Subscript w in (3.8) denotes evaluation of the terms at the mean interface. It is 
noted that = for the Blasius mean flow. 

3.2.2. Wave propagation in anisotropic compliant layers From (2.1), the funda- 
mental wave solutions to (3.1) and (3.2) for a uniformly thick homogeneous layer 
in the transformed variables has the form of 



Tj* = T7* exp[i(a*x* - cot)], 



(3.9) 



where r)* ~ 6cxp(i7a:;| 
and (1.4) yields the following detcrminantal equation: 



For a non-trivial wave, the substitution of (1.7) into (1.3) 



(3.10) 
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where the wavenumbers {k{, k^, k^) = (a*, 0, 7) and 5ps is the Kronecker delta. 



4. Torus translating along an axis of symmetry 

Consider a torus with axes a, 6 (sec figure 2), moving along the z-axis. Symme- 
try considerations imply the following form for the stress function, given in body 
coordinates: 

f{0, V) = (5(V) cos e, (V) sin e, /(V)). (4.1) 

Because of symmetry, one can integrate analytically in the ^-direction obtaining 
a pair of equations for the coefficients /, g in (4.1), 

(sin -ipi — sin ip){a + b cos tp) i 



/(V'l) 



36 

7t(2(a + 6cos V'l))^ Jo (1 — cos(V'i — ^))(2 + a)i 



(4.2) 



2tc 



7r(2(a + ftcosV'i)) 2 Jo 
sin ijji — sin tp 



a + b cos tp 
2 + a 



/(V')[(cos Vi - cos - bl3i)S + PiP] 



cos(V'i — V') 



(sinV^i -sinV>)^ 2 , o 

2 + a - -— - 6 7 5 

1 — cos(V' — tpi) 



+ 



a = a{tp,tp 



(b^ cos V'17 - ^aj F{ln, S)-{2 + a) cos V'i-E(i7t, S) 

62(1 -cos(V'- V-i)) 



(a + 6 cos '(/^) (a + 6cos'!/'i) ' 



/3-/?(^,Vi) = 



dV, (4.3) 

1 — C0S('(/) — t/ji) 

a + b cos ip 

(4.4) 



5. Conclusions 

We have shown how integral equations can be used to solve a particular class of 

problems concerning obstacles in waveguides, namely the Neumann problem for 
bodies symmetric about the centreline of a channel, and two such problems were 
considered in detail. 
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Appendix A. Boundary conditions 

It is convenient for numerical purposes to change the independent variable in (4.1) 

1 1 

to z = y/v^ and to introduce the dependent variable H{z) = {f — y)/v^. Equation 
(4.1) then becomes 

(1 -P)iH + z)H" - (2 + H')H' = H'". (A 1) 
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Boundary conditions to (4.3) follow from (4.2) and the definition of H: 
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if'(0) = -l + eiC7„ + eC4; 




H'{oo) = 0. 



> 



(A 2) 




Appendix B 

A simple sufficient condition for the method of separation of variables to hold for 
the convection problem is derived. This criterion is then shown to be satisfied for 
the ansatz described by (3.27), thus justifying the approach used in § 3. The basic 
ingredient of our argument is contained in the following estimate for a Rayleigh 
Ritz ratio: 

Lemma 1. Let f{z) be a trial function defined on [0, 1]. Let Ai denote the ground- 
state eigenvalue for — d^g/dz^ = Ag, where g must satisfy zLdg/dz + ag = at 
z = 0,1 for some non-negative constant a. Then for any f that is not identically 
zero we have 



Before proving it, we note that the first inequality is the standard variational 
characterization for the eigenvalue Ai . 

Corollary 1. Any non-zero trial function f which satisfies the boundary con- 
dition /(O) = /(I) = always satisfies 
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